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ABSTRACT 

Pulsar timing arrays act to detect gravitational waves by observing the 
small, correlated effect the waves have on pulse arrival times at Earth. This 
effect has conventionally been evaluated assuming the gravitational wave 
phasefronts are planar across the array, an assumption that is valid only 
for sources at distances R 2> 2wL 2 /A, where L is physical extent of the 
array and A the radiation wavelength. In the case of pulsar timing arrays 
(PTAs) the array size is of order the pulsar-Earth distance (kpc) and A is 
of order pc. Correspondingly, for point gravitational wave sources closer 
than ~ 100 Mpc the PTA response is sensitive to the source parallax across 
the pulsar-Earth baseline. Here we evaluate the PTA response to gravita- 
tional wave point sources including the important wavefront curvature ef- 
fects. Taking the wavefront curvature into account the relative amplitude 
and phase of the timing residuals associated with a collection of pulsars 
allows us to measure the distance to, and sky position of, the source. 
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1 INTRODUCTION 



It h as been just ove r thirty years since ISazhinl jl978b 
and lDetweiie3 fl979) showed how passing gravitational 
waves cause disturbances in pulsar pulse arrival times 
that could be used to detect or limit g ravitational wave 
signal strength, and twenty years since iFoster & Backed 
(1990) proposed using correlated timing residuals of a 
collection of pulsars — i.e., a pulsar timing array (PTA) 
— to achieve greater sensitivity. Thought of as a grav- 
itational wave detector a PTA is large: its size L (on 
order kpc) is much greater than the gravitational radia- 
tion wavelength scale A (on order pc) that we use it to 
probe. The conventional approximation of gravitational 
waves as plane-fronted applies only for sources at dis- 
tances R 3> L 2 /X: for less distant sources the curvature 



of the gravitational radiation phasefronts contribute sig- 
nificantly to the detector response. Here we evaluate the 
response of a PTA to discrete gravitational wave sources 
at distances close enough that gravitational wave phase- 
front curvature is important and describe how PTA ob- 
servations of such sources may be used to measure the 
source distance and location on the sky. 

Pulsar timing currently provides the best obser- 
vational evidence for gravitational waves: as of this 
writing the orbital decay of the binary pulsar system 
PSRB1913+16, induced by gravitational wave emission 
and measured by timing measurements, is the most con- 
clusive observational evidence for gravitati onal waves 
dTavlor et alj|l979l : IWeisberg & Tavlonl2005l) . This evi- 
dence is often described, without any intent of dispar- 



© 2010 RAS 



2 X. Deng and L. S. Finn 



agement, as indirect since it is the effect of gravitational 
wave emission — i.e., the effect of the emission on the 
emitter — that is observed^ The effect of the gravita- 
tional radiation, independent of its source, on another 
system — a detecto r — itse l f rema ins to be observed. 

As proposed bv lSazhiJ i 1978) and lDetweilen i ll 9791) 
precision pulsar timing can also be used to detect grav- 
itational waves directly. Gravitational waves passing be- 
tween the pulsar and Earth — i.e., temporal variations 
in the space-time curvature along the pulsar-Earth null 
geodesies — change the time required by successive 
pulses to travel the path from pulsar to Earth. Since 
successive pulses are emitted from the pulsar at regu- 
lar intervals the effect of passing gravitational waves are 
apparent as irregularities in the observed pulse arrival 
times. Through the use of precision timing a collection of 
especially stable pulsars may thus be turned into a galac- 
tic scale gravitational wave detector. 

Timing of individual pulsars has been used to 
bound the gravitational wave power associated with 
stochastic gravitat i onal w aves of cosmolog i cal origin 
teomani & Taylor! 1 19831 : iKaspi et all 1 19941; iLommerJ 
120021) : timing residual correlations among pulsar pairs 
have been used to bound the power in a n isotropic 
stoch astic gravitational wave background jjenet et aU 
2006); the absence of evidence for gravitational waves in 
the timing of a single pulsar has been used to rule-out the 
presence of a proposed supe rmassive black hole binary 
in 3C 66B jjenet et alj|2004l) ; observations from multi- 
ple pulsars have been used to plac e limits on periodic 
gravi tational waves from Sag A* jLommen & Backer! 
l200ll) ; and timing analyses of individual pulsars have 
been combined to place a sky-averaged constraint on 
the merger rate of nearby b lack-hole binaries in the 
early phases of coalescence dYardlev et al.ll201(il) . Ad- 
ditionally, analyses of pulsar timing data have been 
proposed to search for the gravitat i onal wave "mem- 
ory" effect jvan Haasteren & Levinl 120101 : ISetol 120091 ; 
IPshirkov et al.l l201Ch and more ge neral gravitational 
wave bursts jFinn & Lommenl |2010|) . which might be 
generated during, e.g., the periapsis passa ge of a highly 
eccentric binary system (Sesani l2010h or by cos- 



mic string cusps or kinks 1 Damour & Vilenkir] l200ll ; 
ISiemens et al.ll200llLeblond et all2009h . 

Analyses focused on discovering or bounding the 
intensity of a stochastic gravitational wave "background" 
may rightly resolve the radiation into a superposition of 
plane waves with density d s h/dfd 2 Q., where h(t, x) is 
the gravitational wave strain and (/, Q.) denotes a plane 
wave component at frequency / and propagating in di- 
rection Q. Analyses focused on the discovery of gravi- 
tational wave point sources — e.g., supermassive black 
hole binaries — must take proper account of the curva- 
ture of the radiation phase-front, owing to the finite dis- 
tance between the source and the detector, over the detec- 
tor's extent. When detector extent, measured in units of 
the radiation wavelength, is greater than the distance to 
the source, measured in units of detector extent, then cur- 
vature of the radiation wavefront over the detector con- 
tributes significantly to the detector response to the inci- 
dent radiation. 

In [J2]we evaluate the pulse arrive-time disturbance 
associated with a spherically-fronted gravitational wave 
traversing a pulsar-Earth line-of-sight and compare it 
with the response to the same wave in the plane-wave ap- 
proximation. In [J3]we describe how, owing to their sen- 
sitivity to gravitational radiation phasefront curvature, 
PTA observations of point sources can be used to mea- 
sure or place lower bounds on the distance to the source 
distance and, more surprisingly, the distance to the array 
pulsars. Finally we summarize our conclusions in §4\ 



1 Indirect observation is enough to inaugurate the field of gravi- 
tational wave astronomy: i.e., the study of the physical universe 
through the application of the laws and theories to astronomi- 
cal observations involvin g gravitational waves. See, for exam- 
ple. lFinn & Suttoij j2002l) ; ISutton & Find j2002t) , which use the 
effect of gravitational wave emission on several binary pulsar 
systems to bound the mass of the graviton. 
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2 PULSAR TIMING RESIDUALS FROM SPHERICALLY-FRONTED GRAVITATIONAL WAVES 
2.1 Timing residuals 

Focus attention on the electromagnetic field associated with the pulsed emission of a pulsar and denote the field phase, 
at the pulsar, as <po(t). We are interested in the time-dependent phase <po(t) of the electromagnetic field associated with 
the pulsed emission measured at an Earth-based radio telescope, which we write as 

<t>(t) = (j> {t - L - r {t) -TGw(i)] (la) 

where 

(Corrections owing exclusively to the spatial motion of the Earth \ 
within the solar system, the solar system with respect to the pulsar, (lb) 
and electromagnetic wave propagation the interstellar medium J 

tgw = ( Corrections owing exclusively to h(t, X) ) (lc) 

L = (Earth-pulsar distance) . (Id) 

(Note that we work in units where c = G = 1.) 

In the absence of gravitational waves tgw vanishes and the phase front <j>o(t) arrives at Earth at time te(t) = 
t + L + r (t). In the presence of a gravitational wave s i gnal t he phase front arrives at time (t) + tgw (t) ; thus, tgw is 
the gravitational wave timing residual. Following Finn ( 2009) Eqs. (3.26) and (3. 12e), the arrival time correction tgw (t) 
is 

TGw(i) = -^n l n rn "Him(t) (2a) 

where Hi m is the integral of the transverse-traceless metric perturbation over the null geodesic ranging from the pulsar 
to Earth: 



Wlm(t) = L J hi m 



{t + L(„V - L{1 + Z)n) d£, (2b) 



V = (Pulsar location) (2c) 
h = (Unit vector pointing from Earth to pulsar) . (2d) 



2.2 Gravitational waves from a compact source 

At a point X in the perturbative regime far from the source the transverse- traceless gauge gravitational wave metric 
perturbation associated with radiation from a compact source may be written (Finn 1 983) 



h m (t,X)= h lm (f,k)e- 2mft df (3a) 

J —oo 

= f [A lm (f,k)e 2 ^ X - Sl ] e-^df (3b) 

where 

S — (source location) (3c) 
£ 

— tt, FT = (unit vector in direction of wave propagation at X) (3d) 

\X — S\ 

The Fourier coefficient functions Ai m (f, k) are everywhere symmetric, traceless and transverse with respect to X — S: 
i.e., 

A lm (/, fc) = 1+ (/, ~k)e\^ + A* (/, fc)e,W (4) 

where e^' (fc) and e^' (fc) are the usual transverse-traceless gravitational wave polarization tensors for waves traveling 
in direction fc. 
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2.3 Response function 

With expression[3]for him we can evaluate Him in the Fourier domain: 

«h»(/) = L j° h lm (/, V-L(l+ C)n) e - 2 -^ L « 
Aim (f, HO) 



exp [27ri/ + IP - L(l + 0" - <5|)] 



(5a) 



(5b) 



_! \V-L{l + i)n-S\ 

Figure [T] describes the Earth-pulsar-source geometry: note we define R to be the Earth-source distance and 6 the 
angle between the Earth-pulsar and Earth-source lines-of-sight. Along the path traversed by the pulsar's electromagnetic 
pulse phase fronts (V — L£h for — 1 ^ £ ^ 0) we may write 

V-L(l+g)n-S 



k ^ \-p - L(\ 4- - S\ 
r(0 = (distance between gravitational wave phase front and source) 

1/2 



\V- L(l+0n-S\ 



R 



l + 2^f cos 6> 
R 



(6a) 



(6b) 



From Equations [5] and [6b] we note that the time of arrival disturbance owing to a passing gravitational wave de- 
pends on the dimensionless quantities irfL and L/R. PTA pulsar distances L are all on-order kpc. Gravitational waves 
detectable via pulsar timing observations have frequencies greater than the inverse duration of the observational data set 
and less than the typical sampling period: i.e., 10~ 5 Hz > / > 10~ 9 Hz. Strong sources of gravitational waves (e.g., 
supermassive black hole binaries and triplets) in t his band are all expected to be extragalactic, with significant numbers 
within on order 100 Mpc l lSesana & Vecchioll20ld) . The scales of interest are thus either very large or very small: 



e = 10" 
tt/L = 10 4 
ivfLe = 10" 1 



/ 100 Mpc 
lkpc J V R 

/ 



1 yr 1 J \ 1 kpc 



/ 



/ 100 Mpc 
lyr-ijyitepcj V R 

Taking advantage of these scales the Fourier transform Him(f) may be expressed 

r° ~ 

Ui m tf) =L h lm {f,V - L(l + 0") e-' 27TlfLi d£ 



L 



o A lm tf,k(Z) 



i \V-L(l+£)n-S\ 



2nif{-L(+\T-L(l+€)fi-S\) 



(7a) 
(7b) 
(7c) 

(8a) 
(8b) 



= Aim(f, too) 



exp[«/ij(2- i=| 



OS H 
cos 



t/4 



/JR srnO 



erf \e m ^ 4 x 



g (l + COs£>) + l 

1 + cos ' 
fR sin 



1 + cO (c^l. ft 



(8c) 



Except when < (7r/_R) _1/2 < 1 the error function arguments are always large in magnitude. Taking advantage 
of the asymptotic expansion of erf (z) about the point at infinity, 



erf (z) 



exp i 



(9) 
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we may thus write 



where 



and 



A(f,ko) 



Tgw = ~h l h m/ H lm (f) 



(plane wave approximation timing residual) 



= 2A(f, fc )exp 



2ivifL sin 



sine I 27r/Lsin 



(phasefront curvature correction to plane wave timing residual) 

~ r / 2 e 2 

e (1 + cos 9) A(f, ko) exp 7ri/Ll4sin — + — sin ( 

j [n i n m l ;m (/,fc )]exp[27T l / J R] 







nfLe . 2 


): 


sine 


2 Sm 9 \ 



sinc(a;) = 

These expressions are valid for all nfL 3> 1. 



Unnormalized ^ sin(a;) 
sine function 



(10a) 
(10b) 



(10c) 

(lOd) 
(lOe) 

(10f) 
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Figure 1. Gravitational wave source (depicted as a binary system), Earth and pulsar. The gravitational wave phasefronts intersect the 
pulsar-Earth path, traveled by the electromagnetic pulses, in a manner that depends on the distance to the source and the source-Earth- 
pulsar angle 9. 
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3 DISCUSSION 

The Fourier amplitude and phase of the gravitational 
wave timing residual tgw for any particular pulsar de- 
pends on the pulsar's position and distance relative to the 
position and distance of the gravitational wave source. 
Correspondingly, observation in three or more pulsars of 
the gravitational wave timing residuals associated with a 
single source can, in principle, be used to infer the three- 
dimensional source position and location on the sky. In 
this section we discuss the conditions under which suf- 
ficiently accurate measurements of the row are possible 
and provide a proof-of-principle demonstration of how 
they may be combined to determine the sky location and 
distance to a gravitational wave source. 



3.1 Plane wave approximation timing residuals 

Focus attention first on r pw : the plane- wave approxima- 
tion to r gw . In order to relate the angular location of a 
gravitational wave source to the relative amplitude and 
phase of the r pw for two or more PTA pulsars the uncer- 
tainties ol in the Earth-pulsar distances, and <jg in pulsar 
sky locations, must satisfy 



°l < 2pc 



0.1yr _1 \ / 1/2 



/ J V sin 2 (6»/2) 



(Ha) 



»»■> 

The constraint on pulsar angular position accuracy is 
not a severe one; however, the constraint on pulsar 
distance accuracy is quite significant. For this rea- 
son previous pul sar timing searches for periodic grav- 
itational waves dYardlev et alj l20ld; Ijenet et al.1 12004 
lLomme n & Backer 200 ll : lLommer]|200ll) have sought 
only to identify an anomalous periodic contribution to 
pulse arrival times, avoiding explicit use of the ampli- 
tude of r gw . 

Improvements in instrumentation and timing tech- 
niques are making possible pulsar distance measure- 
ments of the accuracy and precision necessary to relate 
the r pw to the gravitational wave source angular loca- 
tion. For example, in recent years a succession of VLBI 
observations have established the distance to J0437 — 
4715 with an accu r acy of ±1.3 p c dVerbiest et alj|2008l ; 
iDeller et alj 120081 ; belled 120091) ; correspondingly, for 
10437 — 4715 predictions of the amplitude and phase 
of Tp W can be made and take part in the analysis of tim- 
ing data to search for gravitational waves at frequen- 
cies / < 0.1 yr. Accurate pulsar distance measurements 



have been identified as a crucial element of a continu- 
ing program to use precision pulsar timing and interfer- 
ome try to test gravity an d measure neutron star proper- 
ties dCordes et alj 120091) ; correspondingly, we can rea- 
sonably expect that the number of pulsars with high- 
precision distances will increase rapidly over the next 
several years. 

Looking toward t he future, observations with the 
SKA dSmits et alj2009h using existing timing techniques 
are expected to be capable of measuring timing paral- 
lax distances to millisecond pulsars a t 20 kpc with an 
accuracy of 20% dTingrav et aifald) . When measured 
by VLBI or timing parallax the fractional uncertainty in 
the pulsar distance is equal to the fractional uncertainty 
in the semi-annual pulse arrival time variation owing to 
the curvature of the pulsar's electromagnetic phasefronts: 
i.e., 



L 



At 



(12a) 



where a T is the timing noise power in the bandwidth 
T -1 (for T the observation duration) about a frequency 
of 2/yr and 



(12b) 



is the peak-to-peak variation in the pulse arrival time 
residual owing to the curvature of the pulsar's electro- 
magnetic phasefronts. Correspondingly, a 20% accuracy 
measurement of the distance to a pulsar at 20 kpc corre- 
sponds to an 0.5% accuracy measurement of the distance 
to a pulsar of the a T at 500 pc: i.e., 2.5 pc. 

This estimate of cjl = 2.5 pc for a pulsar at 500 pc 
assumes constant signal-to-noise and, correspondingly, 
constant timing precision a T . Closer pulsars will have 
greater fluxes and, correspondingly, better timing preci- 
sion; consequently we may regard this as an upper bound 
on the distance measurement error o\l. If we assume that 

is proportional to the pulsar energy flux then 



(13a) 



a T = 0.6 ns ( 

,kpc 



a L = 0.5 pc - — 
\kpc 



(13b) 



i.e., if the only barrier to improved timing precision is 
signal-to-noise then sub-parsec precision distance mea- 
surements should be possible in the SKA era for pulsars 
within a kpc. 

In practice other contributions to the noise bud- 
get will eventually limit the timing precision. Among 
these the most important are intrinsic pulsar timing noise, 
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pulse broadening owing to scattering in the ISM, tro- 
pospheric scintilla tion. Recent im provements in pulsar 
timing techniques dLvne et alJlioiO) suggest that intrin- 
sic pulsar timing noise need not limit timing precision 
measurements in either present or future instrumenta- 
tion. Unresolved pulse broadening owing to scattering 
in the ISM will limit the attainable a T and, thus, ox. 
For nearby (L < 1 kpc) pulsars the measured pulse 
broade ning at v = 1 GHz is typically in the 1-10 ns 
range (Ma nchester et al.ll2005h . Since pulse broadening 
is scales with observation frequency as v~ iA (though 
observations sug gest the index may be closer to —3.9 
dBhatetal.| [2004)) scattering in the ISM also need not 
limit our ability to measure the distance to nearb y pul- 
sars. Finally, recent estimates bv ljenet et al.l d2010 ) indi- 
cate that, for gravitational wave frequencies < yr -1 , tro- 
pospheric scintillation will contribute to the noise budget 
at no more than the 0.03 ns rms level. Consequently, over 
the next decade ox < 1 pc is thus a reasonable expecta- 
tion for millisecond pulsars at distances up to a kpc. 

3.2 Correction to plane wave approximation 

Turn now to the ratio of Fourier amplitudes 

Pif, ko) = 



exp 



irifL ^2 sin 
sin [^nfLes'm 



v e 2 
— I — sin i 
2 2 



(14) 



sin(27r/Lsin 2 §) ' 

When irfLe sin 2 8 ^ 1 each of the sin functions whose 
ratio determines the \p\ has order-unity magnitude and 
\p\ oscillates rapidly from to oo. For a typical PTA pul- 
sar (distance L ~ kpc) this will be the case for source 
frequencies greater than or of order yr -1 and distances R 
less than or of order 100 Mpc (cf. Eqs.|7J. In this regime 
the gravitational wave phasefront curvature plays an es- 
sential role in determining the response of a pulsar timing 
array. 

As an example consider the residuals in the pulse 
arrival times for PSR J1939 + 2134 owing to gravita- 
tional waves from a super massive black hole binary sys- 
tem like that proposed by ISudou et alj d2003h in the ra- 
dio galaxy 3C 66B at RA 2h23mll.4s, De c 42°59'31" 
and l uminosity distance R = 85.8 Mpc (NASA/JPL 
2010). The binary, rule d-out by the subsequent analy- 
sis of ljenetetall j2004), was supposed to have a total 
mass of 5.4 x 10 10 Mq, a mass ratio of 0.1, and a pe- 
riod of 1.05 yr, corresponding to a gravitational wave 



frequency of 1.9 yr . Assuming these binary parame- 
ters to be exact and taking the distance and location of 
J1939 + 2134 as provid ed by the ATNF pulsar catalog 
dManchester et al.ll2005h yields p = 1.1 exp[-0.86«]: 
i.e., the magnitude of the "correction" r cl is greater than 
the "leading-order" plane-wave approximation contribu- 
tion, the magnitude of r gw is 46% of that predicted by 
the plane wave approximation, and the phase of r gw is 
retarded by it/2 radians relative to r pw Q 



3.3 Gravitational wave source distance and 
location on sky 

The magnitude and phase of the plane-wave approxima- 
tion Fourier coefficient function r pw depends on the dis- 
tance to the associated pulsar and the pulsar-Earth-source 
angle. The correction depends on these and, in addition, 
the source distance. Measurements of r gw in a collec- 
tion of pulsars whose angular location and distance are 
known sufficiently accurately can thus be used to mea- 
sure the curvature of the gravitational radiation phase- 
fronts and, thus, the luminosity distance and direction to 
the gravitational wave source. 

Consider a periodic source of gravitational waves 
observed in an array of pulsars whose distances and rel- 
ative locations are known to high precision. The gravi- 
tational wave contribution to the pulse arrival times for 
array pulsar £ may be represented by an amplitude and 
phase, which will depend on the different et = Lt/R, 
nfLe and 6 for each pulsar £ in the array. Requiring that 
these phases and amplitudes all be consistent is a power- 
ful constraint on the source distance and location on the 
celestial sphere. 

As a crude but effective proof-of-principle demon- 
stration of distance measurement by gravitational wave 
timing parallax consider observations made with a selec- 
tion of IPTA pulsars with distances whose precision has 
been projected, as described above, into the SKA era. Fo- 
cus attention on the example source in 3C 66B. For this 
source (7r/) _1 = 0.051 pc. Following Equationll3l four 
pulsars currently monitored as part of the International 
Pulsar Timing Array are close enough that we may ex- 
pect their distance to be measured with a one-sigma un- 
certainty of less than 0.025 pc. Table[Tjlists these pulsars 



2 This comparison is meant only to illustrate the importance of 
r cr in estimating r gw . In this case the distance to J1939+2134 is 
not known accurately enough allow a prediction for an accurate 
estimate of either T pw or t, 



' g w- 
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and their currently measured distances. For these four 
pulsars, 

(i) Calculate the four r gw for the example source in 
3C 66B described above. These complex amplitudes, 
corresponding to periodic timing residual amplitude and 
phase, are our "observations". Denote the r gw for pulsar 

fcbyf fe . 

(ii) Adopt approximate distances to each of these four 
pulsars consistent with normally-distributed measure- 
ment error with standard deviation as given in Equation 

El 



(iii) Using these approximate distances evaluate 



ip 2 (r) 



T k(r) 



(15) 



where r' k (r) is the expected r gw for pulsar k assuming 
the source at distance r and pulsar k at the approximate 
distance found in step \{ii)\ 

The quantity ip 2 (r) is a measure of the misfit between 
the observations Tk and the prediction f' k assuming pul- 
sars at the approximate distances and source at distance 
r. The r that minimizes ip 2 {r) is thus an estimator for 
the distance to the source. Figure [2] shows a set of three 
scatter plots of the estimated r, declination 6, and right 
ascension d> relative to their actual values over 10 4 real- 
izations of pulsar distance errors, for this example. Table 
[2] provides quantitative descriptive statistics for the dis- 
tribution of errors. Clearly, if the pulsar distances can be 
sufficiently accurately determined, the distance and loca- 
tion of a detected periodic gravitational wave source can 
be accurately measured. 

We emphasize that this is no more than a proof-of- 
principle demonstration. It is naive in its treatment of 
uncertainties in the Lk, ignores pulsar timing noise in 
the measured Tk, presumes a sufficiently accurate knowl- 
edge of pulsar locations 9k, and makes use of an ad hoc 
misfit statistic ip' 2 {r). Nevertheless, taken as a proof-of- 
principle demonstration it shows that the potential exists 
for pulsar timing array observations to determine the dis- 
tance to periodic gravitational wave sources if the pulse 
arrival time dependence on the passing waves is properly 
modeled. 



4 CONCLUSION 

Gravitational waves crossing the pulsar-Earth line-of- 
sight lead to a disturbance in the pulsar pulse arrival 





-Z5 -20 
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Figure 2. Scatter plots of the timing parallax estimates of 
the distance, declination and right ascension of a simulated 
gravitational wave source in 3C 66B (i.e., r = 87.8 Mpc, 
dec=2°59'31" and right ascension 2h23mll). See [|33]for de- 
tails. 



time. All previous derivations of this disturbance have as- 
sumed planar gravitational wave phasefronts. The gravi- 
tational wave phasefronts from point gravitational wave 
sources — e.g., coalescing supermassive black hole bi- 
nary systems — are curved, with curvature radius equal 
to the source luminosity distance. The approximation of 
planar wavefronts is thus valid only at distances R 3> 
2nfL 2 /c, where R is the Earth-source distance, L the 
Earth-pulsar distance, and / the gravitational radiation 
frequency. For typical pulsars distances (kpc) and rele- 
vant gravitational wave frequencies (/ < yr _1 ) the cor- 
rection to the disturbance magnitude and phase owing to 
phasefront curvature is thus significant for sources within 
~ 100 Mpc. Here we have derived the curved wave- 
front corrections to the pulsar timing response for such 
"nearby" sources, described when they are important, 
and shown that future gravitational wave observations 
using pulsar timing arrays may be capable of measuring 
luminosity distances to supermassive black hole binary 
systems, and other periodic gravitational wave sources, 
approaching or exceeding 100 Mpc. 

The gravitational wave source distance measure- 
ment described here is properly considered a parallax 
distance measurement. The baselines over which the par- 
allax is measured are, in this case, the timing array 
pulsar-Earth baselines. Crucial to the ability to observe 
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Table 1. Four nearby millisecond pulsars monitored as part of the International Pulsar Timing Array and their current estimated 
distances. These four pulsars are close enough that, in the SKA era, their distances, or the distances to pulsars like them, may become 
known to better than 0.025 pc. 



J-Name kpc 



J0437 


-4715 


0.16 


J1741 


+ 1300 


0.19 


J0030 


+ 0451 


0.30 


J2124 


- 3358 


0.25 



Table 2. Mean, median and standard deviation of the timing parallax estimated distance and location on the sky of a simulated gravi- 
tational wave source in 3C 66B. In this example the source was actually located at r = 87.8 Mpc, dec=2°59'31" and right ascension 
2h23ml Is. Only the four pulsars described in Table[JJwere used to make these estimates. See £|3.3l for details. 



mean 


median 


std. dev. 


r 85.72 Mpc 


85.79 Mpc 


0.99 Mpc 


9 42°59'31" 


42°59'31" 


0°0'5.3" 


2h23m 11.0s 


2h23m 11.0s 


0h0m0.3s 



the effects of gravitational wave phasefront curvature is 
a knowledge of the pulsar-Earth distance L with an un- 
certainty gl < ( 7r /) _1 ~ 2 pc (O.lyr -1 //). We ar- 
gue that pulsar distance measurements of this accuracy, 
while beyond present capabilities (except in the case of 
J0437 — 4715, which is exceptionally bright and close), 
are within the capability of SKA-era observations for 
pulsars at distances L < kpc. 

When considered against the cosmic distance lad- 
der, the distance measurement described here involves 
three rungs to reach distances greater than tens of Mpc: 
first, the determination of the astronomical unit; second, 
the distance to the array pulsars; finally, the distance 
to the gravitational wave source. Since the method de- 
scribed here also provides precise source angular po- 
sition the likelihood of identifying an electromagnetic 
counterpart (e.g., host galaxy) to the gravitational wave 
source is great, raising the possibility of a precise mea- 
surement of both the redshift and luminosity distance to 
a single object at cosmological distances, with the ob- 
vious consequences for independent verification of the 
parameters describing our expanding universe. 

Even in the absence of an electromagnetic counter- 
part measurement of the luminosity distance and angu- 
lar location of a supermassive black hole binary system 

— the most likely source of periodic gravitational waves 

— enables the determination of the system's so-called 
"chirp mass", or (mim2) 3|,5 /(mi + rai) 1 ^ , without 
the need to observe the binary's evolution. 



We have only begun to plumb the potential of grav- 
itational wave observations as a tool of astronomical 
discovery. While this potential will not be realized un- 
til gravitational waves are detected and observations be- 
come, if not routine, more than occasional, explorations 
like these will position us to more readily exploit the op- 
portunities that future observations present. 
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